A Canonical Approach to the Quantization of the Damped Harmonic Oscillator 
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We provide a new canonical approacli for studying tlie quantum meclianical damped liarmonic 
oscillator based on the doubling of degrees of freedom approach. Explicit expressions for La- 
grangians of the elementary modes of the problem, characterising both forward and backward 
time propagations are given. A Hamiltonian analysis, showing the equivalence with the La- 
grangian approach, is also done. Based on this Hamiltonian analysis, the quantization of the 
model is discussed. 



The damped harmonic oscillator (dho) problem is 
characterised by the breaking of time - reversal sym- 
metry. A direct Lagrangian formulation is problem- 
atic because it leads to explicitly time dependent La- 
grangians The standard approach is to 

complement the dho by its time reversed image and 
work with an effective doubled system. The dynami- 
cal group of symmetry of this doubled system is found 
to be SU(1,1) but no unitary irreducible representa- 
tion of the symmetry exists. Time evolution leads 
out of the Hilbert space of states and a satisfactory 
quantization can only be achieved in the framework 
of quantum field theory. This quantization procedure 
is based on the composite Lagrangian of the effective 
system resulting from the doubling of degrees of free- 
dom. The lack of individual Lagrangian prescriptions 
leads to problems in quantization. 

We present in this paper a new method of canoni- 
cal quantization of the dho based on the doubling of 
degrees of freedom. Explicit expressions of the La- 
grangians are given that characterise dual aspects of 
the forward and backward time propagations. We 
have shown that the two cases of overdamped and 
(oscillatory) underdamped motions correspond to dis- 
tinct regimes characterised by real and complex pa- 
rameters, respectively, of the constituent Lagrangians. 
The Hamiltonians corresponding to the complex La- 
grangians are found to be pseudo - hermitian j^. 
We have discussed the diagonalization of the complex 
Hamiltonians pertaining to this regime as a gener- 
alization of the Dirac - Heisenberg method of treat- 
ing the linear harmonic oscillator. The breakdown of 
time - reversal symmetry is manifested in our analysis 
by the appearence of pseudo - hermitian Hamiltoni- 
ans leading to the time evolution of the individual 
modes by nonunitary operators. However, exploiting 



the pseudo - hermiticity of the individual pieces, we 
have shown that well behaved states of the composite 
system are formed. 

We begin with a review of the problem of the 
damped harmonic oscillator (dho). The equation of 
motion of the one - dimensional damped harmonic os- 
cillator is 



mx -\- -\- kx = Q 



(1) 



The parameters m, 7 and k are independent of time. 
If the ratio 



R — — 5- 

4m 



(2) 



is greater than one, the motion is oscillatory with ex- 
ponentially decaying amplitude. Otherwise, the mo- 
tion is nonoscillatory i.e. overdamped. Since the sys- 
tem (|l|) is dissipative a straightforward Lagrangian 
description leading to a consistent canonical quanti- 
zation is not available. To develop a canonical for- 
malism we require to consider (|^) along with its time 
reversed image B 



mij — -\- ky — 



(3) 



so that the composite system is conservative. The 
system (^ and (||) can be derived from the Lagrangian 



mxy -f — {xy — xy) — kxy 



(4) 



where x is the dho coordinate and y corresponds to 
the time - reversed counterpart. Introducing the hy- 
perbolic coordinates xi and X2 § where. 



X = —r^yxi +X2);y = —^{xi - X2) 



V2' 



V2' 



(5) 
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the above Lagrangian can be written in a compact 
notation as 

T ^— ■ ■ -1 ■ 

ij — ^ g^jXiXj i2^ij'^^'^j Q^Qij-^i'^j \y) 

where the pseudo - Eucledian metric gij is given by 
5ii = -522 = 1 and .912 = 0. 

The Lagrangian (^) is invariant under the SU(1,1) 
transformation 



(7) 



where cr is the first Pauh matrix and is an infinites- 
mal parameter. 

The composite Lagrangian (^) is analogous to the 
general bidimensional oscillator Lagrangian 



B . 1^,2 



(8) 



studied recently [g| in connection with the Landau 
problem. Here one exploits dual aspects of the rota- 
tion symmetry of the problem in analysing it in terms 
of opposite chiralities ||]. Symmetry of (|^) under (|^) 
thus offers a possibility of analysing the composite 
theory in terms of systems having opposite chiralities 
w.r.t. the continuous transformations (|^). 

Accordingly we introduce the Lagrangian doublet 



~9ij -^i^j 



(9) 



which are separately invariant under . The Noether 
charges corresponding to the transformations (|^) are 



^ 2 dij^t^j 



(10) 



Thus the systems (O) have opposite "chiralities" w.r.t. 
the transformation (H) which motivates their introduc- 
tion as possible elementary forms of (^. 

The synthesis of Lj^ and L„ is now done by the 
soldering formalism which has found applications in 
various contexts. Duality symmetric electromagnetic 
actions were constructed [|lO| ; implications in higher 
dimensional bosonization were discussed ; the dou- 
blet structure in topologically massive gauge theories 
was revealed fl^ ; a host of phenomena in two dimen- 
sions were analysed However, the analysis that is 
closest in spirit to the one that will be presented here, 
demonstrated the fusion of two one dimensional chiral 
oscillators rotating in opposite directions, into a nor- 
mal two dimensional oscillator 0. Indeed replacing 



gij by 6ij in equation (^) converts it into a doublet of 
chiral oscillators. 



We start from a simple sum 

i(y,z) = L+(y)+i_(z) 

and consider the gauge transformation 

5yi = 5zi = K.i{t) 



(11) 



(12) 



where are some arbitrary functions of time. Under 
these transformations the change in L is given by 



6L{y,z) = 6L+{y)+6L^{z) 

= A,;(j+(y) + Jr(z)) 



where the currents are 



k-^Xi 



(13) 



(14) 



The idea is to iteratively modify L{y, z) by suitably 
introducing auxiliary variables such that the new la- 
grangian is invariant under the transformations (12). 
To this end an auxiliary field Bi transforming as (12), 



5B, = A,; 



(15) 



is introduced and a modified lagrangian is constructed 
as 



L{y, z, B) = i(y, z) - B,{j\+^ [y) + j\ \z)) 
- l{k++k^)B,B, 



(16) 



This lagrangian is now invariant under (p^ ) and (|l5|). 
Since the variable Bi has no independent dynamics, 
it is eliminated by using its equation of motion. The 
residual Lagrangian no longer depends on y or 2 indi- 
vidually but only on the difference y — z. Writing this 
difference as x, the residual Lagrangian reproduces (^ 
with the identification 



r2 



(fc+ + fc-)' 
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k+ + k^ 



(fc+ + fc-) 
(17) 



The essence of the soldering procedure can be un- 
derstood also in the following alternative way. Use 
Xi — yi — Zi in L{y^ z) to eliminate Zi so that 



5u {ViV'j ~ Vt^j - XiVj + ^iXj] (18) 
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Since there is no kinetic term for yi it is really an aux- 
iliary variable. Eliminating yi from L[y,x) by using 
its equation of motion we directly arrive at (|6|) with 
the correspondence ([l7|). Note that the opposite chi- 
rality of the elementary Lagrangians are crucial in the 
cancellation of the time derivative of y in (|l^) which 
in turn is instrumental in the success of the soldering 
method. 

The identification ( |l7| ) has an immediate conse- 
quence. The ratio (0) is found to be, 



„_ fc (fc+ + fc-)^ 



(19) 



For real fc+ , fc_ , the parameters identified by (|l^ cor- 
respond to an overdamped motion of the dhoQ Also 
note that to get the coefficients m and k to be positive 
we require fc_|_ and k- to be of opposite sign, with a 
suitable choice of their absolute values. Finally, for 
positive 7, r > is required. 

Now the physically more important situation is the 
underdamped motion of the dho where the motion is 
oscillatory with decaying amplitude. Here the param- 
eters of (^) are such that the ratio R > 1. As already 
observed this condition cannot be simulated by the 
identification for real values of k±. However, if 
k-^- and k- are continued to complex values so that 



k-\- — K k — — K, 



(20) 



so that the parameters in ( p^ are positive. Substitut- 
ing (|20| ) and (|^ in (^) we get the elementary modes 



T - - i. ■ ^1 

T _ -9 • _ ^ 



(25) 
(26) 



the soldered form of which is the Lagrangian (|^) per- 
taining to the oscillatory limit. Remarkably, the La- 
grangians L± are now complex conjugates of each 
other. 

The Lagrangians (^^ and (^6|) both contain infor- 
mations about forward and backward motions in time. 
To see this we write Lt in the form 



L+ = -igxiX2 - 2 (^1 ~ ^2) 



(27) 



from which the Euler - Lagrange ( E - L ) equations 
follow as 



igX2 — —KXl 

igxi — —KX2 



(28) 
(29) 



According to equations ( p^ , (p^ and ( p^ ) we have 

(30) 



— = n and— = — ^ 
g g 2m 



where 



R=l 



/ Re K 
\ Ini K 



(21) 



then R > 1, which is the required condition for oscil- 
latory motion. Now equation ( [T7| ) gives 



2Re K ' ^ Re k ' 2Re n 

Taking k of the form 

K — Ki + iK2 (23) 

with Ki^2 positive we find that F must be purely imag- 
inary 



(22) 



r = -ig, g>0 



(24) 



*see the discussion below (2) 



TO V 4to 



2 \ \ 2 



(31) 



The solutions to (g^) and ( p9| ) are easy to find. Using 
( pO| ) these solutions can be written in terms of the 
physical parameters of the d.h.o. Now substituting in 
(1), we get. 



X = Aexp (—- — t)exp (iflt) 
2m 

y = Aexp ( — t)exp (— ifJt) 
2m 



(32) 
(33) 



Clearly, x and y correspond to forward and backward 
time propagation with reference to the doubling of 
coordinates (see (0) and @)). The same solutions also 
follow from L_ . In this connection it may be observed 
that Lagrangians structurally similar with (|2^) and 
( p6|) was discussed in Q as the to ^ limit of (||). 
However, it is to be stressed that they are not quite 
identical because the coefficients of (|5|) and (|2^ ) are 
completely different from that of the limiting form of 
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. This is clearly revealed by the calculation of the 
friction coefficient presented in 1^ which comes out to 
be different from that of the actual damped oscillator. 

It will be instructive to look at the problem from 
the Hamiltonian approach. The Hamiltonian follow- 
ing from (^) is 



H = 



1 
2rn 



7 



k 2 W 
2m V 



P2 + '^Xi 



(34) 



where pi = mxi + ^X2, P2 — —mx2 — ^xi are the 
canonical momentum conjugate to xi and X2. Intro- 
duce a canonical transformation from {xi,X2',pi,P2) 
to a;_;p+, where 



/ llj± \i ( milwi \ ^ 



1 



P2 



(35) 



Such transformations, though involving only real pa- 
rameters, were used in Now the composite 
Hamiltonian diagonaliscs as 



H = H. 



where, 



P± 



with the frequencies ll!±, 



2m 



(36) 



(37) 



(38) 



Xi 



=Px ; X2 



the Hamiltonian (B9) becomes 



pI uj\x^ 



(41) 



(42) 



where we have used equations (23), ( {30|) and (3S). 
The above Hamiltonian coincides with of (37). 
Similarly we can derive i7_ from The correspon- 
dence between the Lagrangian and Hamiltonian for- 
mulations is thus complete. 

A question may arise regarding the interpretation 
of the complex Hamiltonians iJ± found in the con- 
stituent pieces. The first point to note is that they 
satisfy 



^1 



(43) 



This hermitian conjugation proprty corresponds to 
the time reversal operation that connects the dou- 
bled degrees of freedom of the closed theory. Also, 
this property manifestly ensures the hermiticity of the 
complete Hamiltonian (^). 

Although iJ± are not hermitian, they are pseudo - 
hermitian 101, 



(44) 



where r\ is the PT operator. Such Hamiltonians have 
occured in the study of PT-symmetric quantum me- 
chanics [Q, in minisuperspace quantum cosmology 
and other constructions |^. To prove the condition 
(H) note that 



r]Xiri 



gijXj 



-gijPj 



(45) 



The Hamiltonians H+ can be shown to follow from 



the Lagrangians L±. Indeed, the Lagrangian (27) is 
already in the first order form . Thus we can read off 
the Hamiltonian directly 



— [x1 
2 ^ ^ 



(39) 



with the symplectic algebra 



The Hamiltonians H±, given by (p7|), are of the 
form 



2 2 

H = ^- + ^-x^ 



(46) 



where x and p are non - hermitian and w is a com- 
plex number. Under rj — PT the operators x and p 
transform as 



\^Xi , Xj } 
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(40) 



From ( |40| ) we find that igxi is canonically conjugate 
to X2- Now by a canonical transformation to the set 
(x,px) defined by 



rjxrj ^ = and r]pr] 

Now define 



uj / ip 



-p 



(47) 



(48) 
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and 



r]\(j)n >= \lpn > 



(59) 



-It "^P 



(49) 



Here a is aptly called the pseudo - hermitian adjoint 
of a with respect to 77. Writing 



we get 



N = aa 



H = uj \ N+- 



(50) 



(51) 



From the basic commutators between the canonical 
variables x and p it is easy to derive that 



Also 



[N, a] = -a 
[N, a] = a 



(52) 



(53) 



Assume that we can construct a complete bidimen- 
sional eigenbasis >, |0n >} diagonalising N 



< (t)n\i^m > = Snm, 



(54) 



Due to (51) this is also the eigenbasis of the Hamilto- 



nian. Now using the commutation relations it can be 
shown that 



N{a\^Pn >) = {n-l)a\ipn > 
Hence we can write 

a\Tpn >= C\i>„-1 > 

where c is some c - number. Similarly 

< 0„|a = d < (t)n-i\ 



(55) 



(56) 



(57) 



The pseudo - hermiticity of N can be exploited to 
relate r]\(f)n > with |-0„ > because 



Nr]\(f>n >= n?7|0„ > 



(58) 



Using the first equation of ( |54| ) we find, upto a phase, 
the following identification. 



The correspondence (|59| ) enables us to reach a crucial 
result 



< </)„|a|'(/)„_i >=< 0„_i|a|-0„ >' 



(60) 

which, along with (|5^) and ( ^ gives 

d = c* (61) 

The last result can be used to show that 

n=\cf (62) 

We find that the eigenvalues of N are real and posi- 
tive. We can also argue that it is integral otherwise 
repeated application of a would yield negative eigen- 
value of N. There, thus, exists a state |0 > which is 
annihilated by a 



a\0 >= 



(63) 



Due to (pll ) this state is the ground state of the Hamil- 
tonian. From the ground state |0 > one can develop 
all the higher energy states by repeated application of 
d. 

From the above solution of the eigenvalue problem 
of (|4^) we can build the physical states of the compos- 
ite system by forming direct product. , Observe that 
due to ( p3| ) the eigenbasis of i7_ will be >, |'0n >} 
if the eigenbasis of is {\ipn >, \<t>n >}• 

Any formulation of the d.h.o. is based on the direct 
[|l| or indirect representation |3|-||]. The direct repre- 
sentation leads to lagrangians having an explicit time 
dependence; hence these are not very popular. The 
indirect representation avoids this problem by a dou- 
bling of the degrees of freedom. It is called indirect 
because, taking the composite Lagrangian and varying 
one degree of freedom yields the equation of motion 
for the other degree (see (^ and its relevant equa- 
tions of motion). The usual composite Lagrangian, 
by construction, is two dimensional. It incorporates 
both forward and backward time propagations. Indi- 
vidual one dimensional Lagrangians displaying these 
properties were non - existant. 

The new point in our paper is that we have pro- 
vided explicit one dimensional lagrangians (equations 
( p5|) and (|2^)) that characterise both forward and 
backward time propagations. Although structurally 
these Lagrangians look two dimensional, the sym- 
plcctic structure effectively reduces one dimension. 
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Moreover we showed that a combination of these La- 
grangians led to (^). In this sense these lagrangians 
are more fundamental. Also they cannot be obtained 
by taking the simple m ^ limit of (||). In the 
region of the parameter space which corresponds to 
the damped oscillatory motion, the parameters of the 
constituent Lagrangians were complex valued. Also, 
these Lagrangians were complex conjugates of one an- 
other. Because of this property, the resulting Hamil- 
tonians were complex valued, satisfying the require- 
ments of pseudohermiticity. This pseudohermiticity 
was exploited to diagonalize the individual Hamilto- 
nians. Based on this, an alternative quantization of 
the damped harmonic oscillator was indicated. 
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